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ABSTRACT 

We present a scenario for non-radiative accretion onto the supermassive black 
hole at the galactic center. Conducting MHD simulations with 1400^ grid zones 
that break the axial and reflection symmetries of earlier investigations and extend 
inward from the Bondi radius, we find a quasi-hydrostatic radial density profile 
p oc r~° ''^ with superadiabatic gradient corresponding to an n ~ 0.72 polytrope. 
Buoyancy generated by magnetic dissipation is resisted by the same fields so ef- 
fectively that energy is advected inward: a state of magnetically-frustrated con- 
vection. This scenario is consistent with observational constraints on energetics 
andouter boundary conditions. 

Subject headings: accretion - magnetohydrodynamics - black hole physics - 
outflows - galaxies: active - methods: numerical 



1. Introduction 

Stellar dynamical measurements indicate a black hole with mass Mbh = 2.4 x IO^Mq 
(Genzel et al. 1997; Chez et al. 1998) at the location of Sgr A* in the galactic center. Apart 
from its origin, a puzzling aspect of this object - and others like it in nearby galaxies - 
is its low X-ray luminosity, given the gaseous environment. High resolution X-ray imaging 
(Baganoff et al. 2001) shows hot gas with temperatures T = 2 keV at densities near = 
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130cm ^ within 1" of Sgr A*. In Bondi (1952) 's theory, gas within the gravitational region 
of influence (Bondi radius) 

_ GMbr , , 

tb = — — ^ 0.03pc (1) 



falls inward, approaching free-fall. This scale is barely resolved by Chandra at the galactic 
center. The natural mass accretion rate is of order AnXr^pCg if the background density is 
P', X — 0.25 for a monatomic gas dominated by thermal (as opposed to magnetic) pressure. 
If matter is converted into radiation at an efficiency rj by the black hole, the luminosity 
L ~ 2 X 10^^(77/10%) erg/s. Reported instead (Baganoff et al. 2001) is a source, potentially 
the hole, with a luminosity of 2.4 x 10^^ erg/s - 10^ times fainter than this estimate. Even 
the advcction of thermal energy across at the Bondi rate would incur ~ 2 x 10^^ crg/s, 
exceeding the observations by three orders of magnitude. Because Bondi assumed spherical 
adiabatic flow of an ideal gas, many factors may be at work in this immense discrepancy. 

First, the flow could be incredibly sporadic and we may have caught it in an off moment. 
However, the dynamical time rs/cg is only ~ 50 years, comparable to the history of radio 
observations and only a few times longer than the X-ray observations. 

Second, it had been suggested (as in the ADAF model of Narayan & Yi 1994) that 
fluid can accrete at the Bondi rate without radiating (77 <C 10%), which might be possible 
if electrons coupled to protons only by Coloumb collisions. However, observations of linear 
submillimetcr polarization (Bower 2000) are interpreted (Mclia ct al. 2000) to imply an 
accretion rate far below Bondi's prediction. Moreover, Bondi-rate influx of electron thermal 
energy alone would exceed hmits by ~ 10^-^. 

Third, Bondi accretion passes through a sonic point only if the effective adiabatic index 
is smaller than 5/3. Gas with 7cff = 5/3 accretes subsonically at all radii, and if 7cff > 5/3 
then quasi- hydrostatic settling flow is expected, becoming sonic just above the Schwarzschild 
radius rsch (Begelman 1978). Monatomic gas with 7 = 5/3 has jeS > 5/3 if viscous or 
magnetic dissipation generates heat (which is not radiated) during inflow (Cowie et al. 1978; 
Scharlemann 1981). Therefore, rapid inflow of the Bondi or ADAF type is unstable to 
motions that reverse the inflow of some fluid elements - potentially reducing Mbh far below 
the Bondi estimate. The ADIOS model of Blandford & Begelman (1999) posits axial outflows 
cancelling the equatorial inflow. Similarly, the CDAF model of Quataert & Gruzinov (2000) 
invokes a rotating convecting atmosphere surrounding the hole. Igumenshchev & Narayan 
(2002, hereafter IN) and Gruzinov (2001) discuss a similar, nonrotating, convecting flow (§3). 

Fourth and fifth, rotation and magnetic fields can both strongly affect the flow and 
interact with one another and with other physical effects in non-trivial ways. 

If the angular frequency is Jls at the Bondi radius, conservation of angular momentum 
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would bring it into orbit at a Kepler radius rx — figr^/Cj. If cooling permits the forma- 
tion of a thin disk, this signals the onset of efficient radiation (7/ ~ 10% as used above). 
If angular momentum transport is weak, however, this can be a limiting step in accretion. 
Nayakshin (2003) has recently suggested a cold disk onto which the hot flow condenses, 
accreting sporadically with a long duty cycle. The disk of stars inferred by Levin & Be- 
loborodov (2003) could have arisen this way. Recall, however, condensation onto this disk 
would exceed constraints if it emitted X-rays. 

Similarly, magnetic fields can grow in strength as they arc dragged inward by the flow. 
Accretion will shear embedded flelds to become radial. An initial fleld, whether uniform or 
tangled (Shapiro 1973), will be stretched toward the split-monopolc (hedgehog) configuration 
in which oc (e.g., IN), incurring a centrally divergent energy. Although radial fields 
exert no net force, this configuration is unstable and unattainable (Zeldovich & Novikov 
1971). Instead, infiow can stall at a magnetic turnaround radius Tmag 
where VA{r) is the Alfven velocity. Mass inflow will be limited by the rate of magnetic 
energy dissipation via reconnection, by interchange instabilities, or by infiow along open 
field lines. 

This contrasts with the widely-held picture of magnetic fields enhancing accretion 
through angular momentum transport, either via turbulent stresses (Balbus & Hawley 1991) 
or via magnetoccntrifugal winds (Blandford & Payne 1982). Also, if fields are generated lo- 
cally through convection or magnctorotational instability, then they are limited in strength 
to partial cquipartition with kinetic energy. Fields strengthened by inflow can reach equipar- 
tition with the gravitational potential, which can be a larger value (§2.1). 

We pause here to note the critical importance of the density profile in the phenomenology 
and viability of models; see also Quataert (2003) and Gruzinov (2001). An atmosphere with 
p oc r~" and, generically, T oc r^^, has a bolometric free-free luminosity L oc r^p^T^/'^ oc 
j,(5-4n)/2 ^Qj, ^ single frequency, Lj, oc r^p^T^^^"^ oc 7''''/2-2n^^ long as T > 3 keV. Thus 
profiles with n > 5/4 (or n > 7/4) suffer a singularity in the bolometric (or single- frequency) 
luminosity. If n < 5/4, emission is dominated at outer radii; this is consistent with the small 
contrast in Chandra images. 

Similarly, the thermal time varies as T^/'^ / p oc r""-*^/^ whereas the infiow time r/v — 
Anr^p/MBB. oc r^~". Thus, fiows with n < 7/4 are dynamically nonradiative. Finally, 
the local Bondi accretion rate varies as r'^pCg oc r^/^~": fiows with n < 3/2 accrete at 
M ~ MB(i?B/i?Sch)=^/'-". 

Hydrostatic models are characterized by a polytropic index n — l/(7efr — 1), which 
describes the correlation between density and pressure p oc p'^^^ — In a Keplerian 
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potential, p oc r"" as used above. When — 5/3, both hydrostatic and Bondi-hke (e.g., 
ADAF) flows have n = 3/2, shallow enough to be advective but too steep to avoid the 
luminosity constraint. To force a shallower density proflle requires additional pressure, either 
through rotation, magnetic fields, or an entropy inversion. The latter is convectively unstable. 

CDAFs pass this test by achieving n = 1/2 in a combination of rotational, turbulent, 
and thermal support. However, the Quatacrt & Gruzinov (2000) model requires the inward 
angular momentum transport due to buoyant convection to exceed the outward transport 
due to magnetic fields. This may be possible for fields due solely to magnetorotational 
instability (Narayan et al. 2002), but it is unlikely a property of field strengthened by shear 
in infiow. 

More seriously, any model invoking rotational support is plagued by outer boundary 
conditions. Outside one expects solid-body rotation on radial shells, implying low specific 
angular momentum (j) near the axis. Sufficiently low-j gas can fall directly within the 
innermost stable orbit; this comprises a fraction ~ (r^/rx)^'^^ of the total. If j is somewhat 
higher, fluid falls inward to its own Kepler radius and shocks; the resulting pressure gradients 
drive a quadrupolar outflow of low-j material along the equator. This has been observed 
in nonmagnetic simulations by Proga & Begelman (2003) and by us. It could be avoided if 
the inner boundary supplied a strong jet that interfered with axial accretion, as occurs in 
protostars (Matzner & McKee 2000), but suggestions of such behavior in Chandra images are 
recent and tentative (Morris ct al. 2002). Similarly, models that invoke rotational support 
to stem accretion require that rotational support remains important all the way out to Tb - 
implying an asymmetry to the X-ray images which is not observed. 

How can gas establish the flattened density proflle required for low-luminosity accretion, 
while being fed low-angular-momcntum, lightly magnetized {j3 > 10 for 5 < 1 mG) material 
at r^? IN and Gruzinov (2001) offer a solution to this conundrum (§3); our simulations 
suggest a different answer. 

2. Simulations 

In order to further explore the interaction between infall, rotation, magnetic flelds, and 
buoyancy, we conduct MHD simulations with 1400'^ zones arrayed in a uniform Cartesian 
grid, the largest MHD simulations to date. These were performed on the CITA McKenzie 
cluster: 512 Pentium-4 Xeon processors running at 2.4 GHz (Dubinski et al. 2003). At 
this resolution, each full dimensional sweep corresponding to two timesteps took 40 seconds. 
The code (Pen et al. 2003) is based on a 2nd order accurate (in space and time) high 
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resolution Total- Variation-Diminishing (TVD) algorithm. It exphcitly conserves the sum of 
kinetic, thermal and magnetic energy; hence magnetic dissipation (at the grid scale) heats 
gas directly. No explicit resistivity or viscosity is added, and reconnection and shocks occur 
through the solution of the flux conservation laws and the TVD constraints. Magnetic flux 
is conserved to machine precision by storing fluxes perpendicular to each cell face. 

Inner boundary conditions are imposed on a cube of width 24 grid cells. Interior to 
the largest inscribed sphere within this cube, gravitational forces are turned off. At the end 
of each time step magnetic flelds in the cubical region are relaxed to the vacuum solution, 
permitting rapid reconnection in the interior zone. The Alfvcn speed is matched to the 
circular speed at the surface of this region by the removal of matter, and the sound speed is 
matched to the same value by the adjustment of temperature. The pressure of this matter 
is always smaller than that of infalling material, and we never observe spurious outflows. 

Whereas many prior simulations have started with a disklike, rotation dominated ge- 
ometry and no \ow-j material, we wished to investigate the potential role of axial infall. We 
attempt to separate generic physical effects from artefacts of boundary conditions. This is 
challenging, since the energy available at the center of the simulation is always larger than 
in any other region. Results may depend sensitively on the choice of inner boundary, and it 
is not feasible to simulate directly from to rsch- We search instead for scaling relations 
connecting these radii, bearing in mind the possibility of strong reactions such as outffows 
from an unresolved inner region. 

For this reason, and since we are interested primarily in the effect of external parameters, 
most of our runs have worked inward from the Bondi radius. Simulations far inside tb 
elucidate local phenomena, but cannot easily be matched to their exterior environment. At 
each step, the outer 20 grid cells are replaced with the values from the initial conditions 
to enforce the continued inflow of new material. We avoid threading the outer boundary 
with any magnetic flux by adding flux lines to the fluid inside a region 1 / 2-3 /4 of the box 
size, with an admixture of random field loops and large loops that thread the whole box. 
(Most of the energy is in the coherent loops.) We maintain a mean value /3 = 10 in the fiux 
generation region. To estimate Tmag; we assume B oc as in the split-monopole. This 
gives Tmag — 100 grid cells. The estimate is not rigorous: fields are added in a cubical region; 
rotation changes the shearing rate, and B oc is not perfectly achieved. 

To avoid immediate energetic feedback from the central regions due to the initial condi- 
tions, we start with an empty interior (an evacuated sphere of half the box radius), and let 
matter fall in. We then watch the evolution away from Bondi to draw conclusions about the 
fate of the flow. In grid units, box length is 1400, and GM = 700. We set Cs{r = oo) = 1, 
and p{r = oo) = 1. The Bondi radius was set at = 700, touching the closest box edge. 
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The simulations ran successfully for 6000 time steps to t — 650, or 1.5 free fall times from 
vb- At that point, magnetic fields squeezed out along the midplane to the outer boundary, 
leading to numerical instabilities. 

With a resolution of 10^-^^ radial zones, of which 10^-^^ is used for central and outer 
boundary conditions, we have 1.34 decades of scale in which to arrange the Kepler radius 
rx and the magnetic turnaround radius rmag- We have performed runs that arc initially 
rotationally supported (rmag < "f^K) and purely rotationally supported ones (r^ag = 0), as 
well as magnetized but nonrotating initial conditions {r^ = 0). The fluid is initialized with 
solid body rotation on shells, and constant specific angular momentum on cones of constant 
polar angle. The production simulation (fig. 1) chose rx = 336, rmag = 100. To further 
break any discrete symmetries, the velocity field was modulated at the 5% level at multipoles 
up to Z = 2. 

These simulations differ from previous ones in two important ways. First, magnetic and 
viscous dissipation occur only at the grid scale (or inner boundary); no enhanced diffusivity 
or viscosity is applied. Second, we break the alignment usually assumed between magnetic 
and rotational symmetry axes. This was accomplished by introducing large-scale flux loops, 
misaligned with respect to the rotational axis, that are dragged toward the central object. 

2.1. Results 

Simulations were run about one dynamical time at rg. Although this does not give us 
any handle on the long-term evolution of the flow once its inner behavior begins to alter the 
conditions outside r^, it does yield a picture of how the inner flow responds - for many inner 
dynamical times - to initial and boundary conditions. A snapshot of the resulting magnetic 
field structure is shown in 1, and various stresses balancing gravity are plotted versus radius 
in figure 2. 




Fig. 1. — Magnetic field structure: midplane (x — y) slice. Magnetic pressure is shown along 
with projected magnetic field vectors. The inner half of this picture corresponds to figure 2. 
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Fig. 2. — Run of supporting stresses (radial components of momentum fluxes), averaged on 
radial shells and compared to the local gravitational force per unit volume. Also plotted is 
the standard deviation of gravitational force per unit volume on radial shells. 
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At this point in the simulation, a central hydrostatic region supported by gas pressure 
- not rotation or magnetic fields - has grown outward. Within the region plotted in figure 
2, p oc r"" for n ~ 0.72 also, P oc r~^'^^ (so that T oc r~°'^^, whereas was expected), and 
magnetic pressure Pmag ~ Gas pressure gradients dominate magnetic stresses by a 

factor ~ 10; Reynolds stresses, including rotation and inflow, are smaller still (and switch 
direction). Rotation is roughly one tenth the Kepler rate, despite that the entire region is a 
factor of two inside the initial tk- The ratio of Alfven to inflow velocities is similarly ~ 10, 
indicative of magnetic braking. 

In contrast to Bondi and ADAF-type flows, inflow is very subsonic; correspondingly, 
the mass accretion rate is strongly suppressed relative to Bondi's estimate. However, the 
accretion rate does agree with the Bondi rate derived from conditions at the inner bound- 
ary. This state is not rotationally-supported like the CDAF. It resembles in some ways the 
"CDBF" of IN and Gruzinov (2001); but see §3 for differences. 

To test the role of magnetic fields in the clogged infiow, we suddenly turned the magnetic 
fields off, and evolved the fiuid for a dynamical time. The fiow returned to Bondi's solution. 



3. Physical Interpretation 

The fluid can be modelled as one in quasi-hydrostatic equilibrium with a polytropic 
index 7eff ~ 2.25. Since the adiabatic index is 7 = 5/3, this represents a strongly supera- 
diabatic state. In the usual description of entropy-driven convection, this can only occur 
when convective velocities approach the sound speed, which, in a power-law atmosphere, is 
roughly free-fall. 

Saturation at a constant Mach number is a feature of Quataert & Gruzinov (2000) 's 
CDAF model and the CDBF model of IN, both of which have n — 1/2 and fall within 
class II of Gruzinov (2001) 's classification of self-similar nonradiative fiows. This slope 
is clearly fiat enough to satisfy all of the observational constraints discussed in §1. The 
value of n = 1/2 derives from assuming a positive convective luminsoity Lconv ~ i-e., that 
the gravitational energy released bubbles out through the fiow rather than being dragged 
inward (as it is in Bondi and ADAF flows). If Cg scales with the Kepler speed Vk-, then 
Lcorw r'^pv\AA oc r^/^^"A^, where M. is the turbulent Mach number. Constancy of Lconv, 
required by steady state, implies M. oc r"^~^/^. If n > 1/2, rises outward until it saturates 
at ~ 1; thereafter, n = 1/2. 

Our simulation shows n = 0.72 and mildly subsonic convective motions: A4 ~ 0.4 - 
not increasing outward at r°'^^ as implied by the above argument. This has prompted us to 
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calculate the total flow (Lconv) of magnetic, thermal, and gravitational energy. We flnd it to 
be inward a,n.d small: -Lconv — —0.04 x (Anr'^pvj^). 

How can subsonic convection persist in the presence of strong superadiabatic gradi- 
ents? And, how can the net flow of energy be inward despite the existence of convective 
motions? Since the detailed treatment of magnetic fleld structures is the only physical ef- 
fect distinguishing our final state from ordinary gaseous convection, the simulation results 
must represent a state of magnetically-frustrated convection in which magnetic shear stresses 
oppose buoyant motions. This is verified by the cquipartition between magnetic stresses 
and buoyant stresses (variations of pg at fixed r) seen in figure 2. Correspondingly, kinetic 
energy is suppressed by a factor of three compared to the free energy available in density 
fluctuations (variations of pGM/r at flxed r). For further corroboration we note an anticor- 
relation between magnetic and buoyant stresses (cross-correlation coefficient -0.2) and also 
the relaxation to Bondi flow when flelds were removed. An analogy to magnetic frustration 
is low-Reynolds-numbcr convection, in which viscosity slows or halts buoyancy. (That said, 
energy inflow may be a feature of our inner boundary conditions - a question for further 
study.) 

The coincidence in slope between gas and magnetic pressures most likely arises from 
the growth of magnetic flelds to balance buoyancy: since buoyancy scales with gas pressure 
in the presence of a strong superadiabatic gradient, so must the magnetic stress. 

One can discuss the plausible scenarios of the global gasdynamics of the Sgr A* region. 
In addition to the ambient gas, stellar winds may inject a larger mass flux into the region 
than would be lost by even ideal Bondi accretion (Quataert et al. 1999). If the natural 

accretion rate is sufficiently slow, the injected mass could conceivably push the ambient gas 
outwards. Similarly, the cooling time at the Bondi radius is approximately 10^ years, and 
on this time scale matter has to either flow in, or be pushed out. 

4. Conclusions 

We have presented a new picture of the hot plasma in the vicinity of the galactic cen- 
ter black hole based on the results of new very large three dimensional MHD simulations. 
We find that the code produces results consistent with observational constraints, including 
boundary conditions and total luminosity. Fluid remains in quasi-static equilibrium sup- 
ported primarily by thermal pressure, with a radial density profile p oc r~^^'^'^ in reasonably 
spherical symmetry. This implies a significant entropy inversion; we believe buoyancy is 
impeded by subdominant magnetic fields. 
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A distinguishing feature of the flow is that energy is advected inward, despite convection, 
albeit quite slowly compared to the local estimate 47rr^p(r)f /^(r)^. The same is true of central 
regions of Bondi flow, where E = —Mcs{r > / (7 ~ 1)- Unlike Bondi, the M is also much 
slower than the dynamical accretion rate. 

This feature provides a point of contrast with related suggestions for the Sgr A* flow: 
CDAF (Quataert & Gruzinov 2000) and CDBF (IN), both of which obey p oc r'^''^ and 
contain different levels of rotational support. The 1/2-law in these models is a consequence 
of a positive outward convective luminosity, plus the saturation of convective motions at 
constant Mach number (sec § 3 and Gruzinov 2001). However, the density profile is essentially 
unconstrained when E and M (as well as the angular momentum flux J) are all effectively 
zero (Gruzinov 2001, his class- IV flow). 

The influx of energy also provides a method of distinguishing magnetically-frustrated 
flows from the other models even if the central density proflle cannot be probed. In CDAFs 
and CDBFs, energy must be convected outward and will heat gas outside the Bondi radius 
at a rate comparable to Mc^ ~ 1O^L0. If present, effects of this heating on the hot gas 
surrounding the hole may be visible. On the other hand, winds from massive stars in the 
same volume may overwhelm inward advection. 

Blandford & Begelman (1999) criticize ADAFs on the basis that they possess a pos- 
itive BernouUi function [B > 0), and that this predisposes them to outffow (reahzed in 
their ADIOS model). Note however that B is positive (though small) in Bondi ffow, where 
it is set by the external pressure; however, no outffow develops. We find ;B > in the 
magnetically-frustrated ffow as well. We do observe upwelling, but it takes the form of an 
overall quadrupolar circulation with inffow along the axis and outffow along the equator. 
There is no indication that this upwelling becomes supersonic at any radius. 

We speculate that non-radiative MHD ffows are in general very inefficient at accreting. 
The corollary is that accretion only occurs when the cooling time is short. 

This research was funded by NSERC, by the Canadian Foundation for Innovation, and 
by the Canada Research Chairs program (for CDM). We thank Phil Arras, David Ballantyne, 
Robert Fisher, Chris McKee, and an anonymous referee for comments. 
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